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Our procedure to obtain best values for cosmological parameters from five recent multiparameter 
fits is as follows. We first study the values quoted for r, a s , lu + l(wo + l), wi and flfc, arriving at the 
conclusion that they do not differ significantly from zero, and their correlations to other parameters 
are insignificant. In what follows they can be therefore ignored. The neutrino mass sum T,m v also 
does not differ significantly from zero, but since neutrinos are massive their sum must be included 
as a free parameter. We then compare the values obtained in five large flat-space determinations of 
the parameters £m„, uit, Mm, h, r, n s , A s and ag. For these we compute the medians and the 
f 7-percentile and 83-percentile errors by a described procedure. 



PACS numbers: 98.80.-k, 98.80.Es, 95.35.+d 

I. INTRODUCTION 

Since July 2004 five large analyses HUH HQ of cos- 
mological parameters have appeared (we call them pa- 
pers PI, P2, P3, P4 and P5) based on partly overlapping 
data. They represent a wealth of new information which, 
however, may appear contradictory to the outsider and 
difficult to digest. Since each analysis determines param- 
eter values separately for different subsets of data and for 
several choices of parameter spaces, the results are legio. 
For instance, these five papers quote altogether 56 differ- 
ent values for the density of matter in the Universe, w m . 
The parameter values exhibit notable differences due to 
differing choices of priors, methods of analysis, simula- 
tion variance and statistical fluctuations. The purpose 
of this review is to compare these five analyses, and to 
try to extract as reliable values as possible. Although 
the authors of the analyses may be unlikely to believe in 
the results of anyone else, theoreticians do need recom- 
mended values. 



II. DATA SETS 

All analyses quote results from fits to different sub- 
sets and combinations of data sets. This is of course 
of importance in order to test different sets for consis- 
tency or for correlations. However, since the purpose of 
the present analysis is to obtain as accurate information 
about parameter values as possible, we consider only fits 
to maximal data sets. 

All the five papers P1-P5 use the CMB (TT) Q and 
(TE) power spectra on large angular scales from 
the first year Wilkinson Microwave Anisotropy Probe 
(WMAP) observations. PI .lj combines this with the 
constraints from the Sloan Digtal Sky Survey (SDSS) 
galaxy clustering analysis [8j , the SDSS galax y b ias anal- 
ysis at z = 0.1 and the SNla constraints [Hj]. A dis- 
tinctive feature of the PI analysis is the inclusion of Ly-a 
forest data [Tlj . PI mentions explicitly that inclusion of 
power spectra from the CMB observations on small an- 
gular scales by the Cosmic Background Imager (CBI) in 



2004 [13, the Very Small Array (VSA) pj, and the Arc- 
minute Cosmology Bolometer Array Receiver (ACBAR) 
[fH do not affect their results much, nor does the inclu- 
sion of the 2 degree Field Galaxy Redshift Survey (2dF- 
GRS) 15j galaxy clustering power spectrum. 

P2 uses an SDSS "Main sample" H3 of galaxy 
clustering data which is of earlier date than the sam- 
ple used in PI. The distinctive feature of the P2 analysis 
is the inclusion of spectroscopic data for a recent sam- 
ple of 46 748 luminous red galaxies (LRG) in the range 
0.16 < z < 0.47 which show a significant acoustic peak in 
the redshift-space two-point correlation function (related 
to the galaxy power spectrum). P2-P5 do not include 
Ly-a forest data, some of them mentioning that the in- 
formation is still susceptible to systematic errors. 

P3 3 uses in addition to WMAP d the follow- 
ing sets of CMB data: the BOOMERANG (TT) power 
spectra, the BOOMERANG (TE) + (EE) + (BB) po- 
larization spectra, and the (TT) power spectra of DASI 
[13, VSA mi, ACBAR []j, MAXIMA [3, and CBI 
2004 data O. The shapes of the SDSS and 2dFGRS 
"T, 1 1 olj are used as constraints, but not 
In addition the SNla "gold" set [13 is 



power spectra 
the amplitudes 
included. 



P4 [4| relies on either the 2dFGRS [15J or the SDSS [8j 
ower spectrum of galaxy clustering , the WMAP data 
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but in contrast to PI they also use the (TT) power 
sjjectra of VSA [13, ACBAR 0, and CBI 2004 data 
UJ. We use the parameter values derived from 2dFGRS 
data rather than SDSS, because P4 represents 2dFGRS. 

P5JJ analyzes exclusively CMB data: (TT) from 
m l llT j and in particular all the(TE) + (EE) polarization 
data now available 

II Hlllla They do not use LSS 

constraints. 

It would have been interesting to include also the re- 
cent analysis of cosmic shear [jjj. However this would 
not have been consistent because the parameter space 
analyzed is quite different from that of papers PI - P5. 
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III. PARAMETER SPACES 

The most general parameter space explored in the five 
papers is 13-dimensional 

p = (uj b ,u) m ,VL kl Y>m u ,T, h, a s ,b,w,Wi,n s ,A a ,a s ,r) . 

(1) 

The normalized matter density parameter f2 m does not 
appear explicitly in this list because it is given by Ct m = 
u) m /h 2 where h is the Hubble parameter in units of 100 
km s _1 Mpc -1 . Similarly, fi;, = u>b/h 2 is the normal- 
ized baryon density parameter. One has u> m — ujdm ± ^>b 
and hence Q m — ^dm + ^f>, where fW = uJdm/h 2 — 
Qcdm + Qv is the density parameter of dark matter, some- 
times denoted J7 C for CDM. The density parameter of 
dark energy is then Q\ = 1 — Q m — Qfc, where fife is the 
curvature or vacuum density. The ratio of pressure to 
energy density for dark energy is w = wo + it>i(l — a), 
where a is the scale parameter of the Universe, erg is the 
rms linear mass perturbation in 8/i _1 Mpc spheres. 

The fraction of dark matter that amounts to neutrinos 
is fv = Qv/Qdrn- Assuming that the neutrinos are Ma- 
jorana particles with standard freeze-out, the sum of the 
three neutrino masses is £m„ = 94.4 u>dmfv eV. 

The parameters r, n s , a s , A s , r, b describe fluctua- 
tion properties: the scalar spectral index is n s , the scalar 
amplitude is A s , the tensor spectral index n t is not in- 
cluded because there is no information, the ratio of tensor 
to scalar amplitudes is r = A t /A s , and the running of n s 
with k is a s . The parameter r measures the Thomson 
scattering optical depth to decoupling, and b is the bias 
factor describing the difference in amplitude between the 
galaxy power spectrum and that of the underlying dark 
matter. For b no value is quoted. 



IV. METHOD OF ANALYSIS 

Since the data sets are to a large extent overlapping 
they are not independent. Moreover, parts of the re- 
ported errors are systematic and not stochastic, so dif- 
ferent data sets cannot be summarized using statistical 
methods, either frequentist or Bayesian. The breakdown 
into data subsets differing in methods of data sampling, 
statistical analysis, and simulation demonstrate notable 
variations in the parameter values (see for example Chu 
et al. 20)). We clearly want to use results from as large 
data compilations as possible, not from subsets. 

Many parameters are determined with quite low pre- 
cision, the published la errors are large, the correlations 
likewise insignificant, and the dependences on the choices 
of data subsets, priors, and parameter spaces are statisti- 
cally weak. The tendency in the cosmological literature is 
to draw optimistic conclusions from la and 2a confidence 
regions obtained by marginalization in many-parameter 
spaces. We think that parameter estimates from the five 
data sets, at their present state of accuracy, can well be 



TABLE I: Comparison of empirically defined errors 



a± ± Sa± 


ai ± 5a,2 


Our method 


\a 2 - ai\/2 
+5ai 


1±5.0 


2 ±5.0 


1.5 ±4.80 


1.5 ±5.5 


1 ± 1.0 


2 ±1.0 


1.5 ± 1.08 


1.5 ± 1.5 


1±0.5 


2 ±0.5 


1.5 ±0.71 


1.5 ± 1.0 


1±0.2 


2 ±0.2 


1.5 ±0.58 


1.5 ±0.7 


1±0.5 


2 ±1.0 


-, on +1.09 
l.OO _ .67 





combined without sophisticated methods. We then pro- 
ceed as follows. 

Each measurement i of a parameter a we represent by 
a Gaussian ideogram with central value ai, symmetrized 
error Sai , and unit area [in least squares fitting one weighs 
the area by (Sai)~ 2 ]. We sum these ideograms, and de- 
fine our best values as medians (as do for instance PI 
and P5), not means. In lieu of variances we define our 
errors as the 17-percentiles and the 83-percentiles of the 
sum ideogram. The most notable effect is then that these 
" errors" are always larger than the smallest individual er- 
ror in the data set, whereas in statistical averaging errors 
always come out smaller than the smallest error entering. 
One may take this as an estimate of systematic errors for 
which there anyway does not exist any statistical pre- 
scription. 

Some numerical examples in Table I demonstrate the 
properties of this method. In comparison we give the re- 
sults for another empirical method in which |o,2 — a\\/2 
is taken as an estimate of a systematic error which is 
added linearly to the statistical error 5a\ = 5a2- (An- 
other prescription is needed for asymmetric errors.) The 
last line in the Table exemplifies the combination of data 
with unequal errors by our method. The penultimate line 
exemplifies data which clearly are in contradiction. No 
such flagrant contradictions appear in the data that we 
shall combine. 

Of course no probability measure or statement of sta- 
tistical confidence can be associated with our type of er- 
rors. Whenever wc mention la errors or confidence limits 
(CL), they refer to quantities in the data sets. 

V. PARAMETER VALUES 

A. r, a a , w, w , Wi, Q k 

Let us first discuss parameters for which the theoretical 
motivation is model dependent and the evidence insignif- 
icant. 

PI quotes three upper limits (95% CL) for the parame- 
ter r, the smallest value being r < 0.36 in a 7-parameter 
fit with a s — 0, £to„ = and w = —1. The effect 
of including r as a seventh free parameter is that the 
information on the basic six parameters is degraded - 
their errors increase. Their central values change only 



3 



within the quoted la errors, thus r is not significantly 
correlated to the basic parameters. The largest correla- 
tion is with n s , still an increase of only about la (this 
is discussed in detail in PI p|). When any of the pa- 
rameters a Sl Em„, w is allowed to vary, the parameter 
space becomes 8- or 9-dimensional resulting in a further 
deterioration of the information on the basic parameters 
as well as on r. 

P3 quotes r < 0.36 and P4 quotes r < 0.41 from fits 
in the same 7-parameter space. The inclusion of r as a 
free parameter degrades slightly the information on some 
basic parameters. As in the case of PI, r is not signifi- 
cantly correlated to any other parameter except slightly 
to n s , at the level of about la. 

P2 and P5 do not quote any value for r. Thus there 
is no significant evidence for r ^ 0, only for a commonly 
acceptable limit r < 0.4. 

The parameter a s is determined by PI and P3. PI 
quotes fits in two different parameter spaces: 



r < 0.45, a s 



-0.006 



+0.012 
-0.011' 



T, and 



I +0.08 



r < 0.45, a s = -0.011 ± 0.012, w = -0.91l£ 09 . 
In neither case is a s significantly different from zero. 

P3 quotes a s — — 0.051 j^^g when w = — 1 in rather 
marked disagreement with PI. Because of this conflict 
we choose not to combine it with PI. The effect of in- 
cluding a s in the PI and P3 fits increases the errors of 
the basic parameters considerably. In PI no significant 
correlations are found, in P3 there are some correlations 
at the la level. That the running of n s is poorly deter- 
mined is not surprising since we shall see later that also 
n s is poorly determined. 

To proceed with w we take r = and a s = 0. With 
this assumption four values of w have been determined: 
PI: -0.99 ±0.09 , P2: -0.80 ±0.18 , P3: -0.94±°:io. 
P4: — 0.85±g;if. 

The inclusion of w as a seventh parameter in PI and 
P3 has very little effect on the errors and central values of 
the basic six, more so in P4, but the correlations are still 
insignificant. One notices in PI, however, that keeping r 
and a s free improves the w errors and reduces its central 
value by slightly more than la. 

Including w as a fifth parameter in P2, degrades the 
errors of their basic four parameters, and changes their 
central values by less than la. 

The above values can then be used to define the median 
and errors 



w = -0.93 



+0.13 
-0.10 



(2) 



PI also determines the parameters wq, w\ in the com- 
bination w = Wq + it>i(l — a), with the result 
w = -0.98 ± 0.19, wi = Omto H- 
Thus there is no significant information indicating w ^ 
— 1 or w\ ^ at present. In a Bayesian study of the need 
for either of the parameters w or n s (using a modified 
PI data set) Mukherjee et al. conclude that theories 
with more parameters are not disfavored, but that the 
improvement does not warrant the additional complexity 



in the theory. Thus we choose results from 7-parameter 
fits which assume w = — 1. A fortiori w\ is useless. 

The available determinations of the curvature para- 
meter Ofc for the case w = — 1 are 

P2: -0.010 ± 0.009, P3: -0.027 ± 0.016, P4: 

-0 74+ 049 
u.ud_ 052 . 

The general tendency is confirmed that the inclusion 
of flk degrades the information on the basic parameters, 
also the correlations are insignificant with one exception: 
there is a positive correlation between fi& and h, most 
visible in P3. The above set of data define the median 
and errors 



fin 



i-n 



-0 023+ 017 

U.UZO_g 050 . 



(3) 



We conclude that there is very little information on 
curvature at present, and that h to some extent takes its 
role. 

In the sequel we shall only use parameter values deter- 
mined under the assumptions of flat space and r = a s = 
w + 1 = 0. This maximizes the information on the basic 
parameters, and it does not entirely neglect the effects of 
possibly non-vanishing fi^ or r because of the presence 
of the slight (fife, h) and (r, n s ) correlations. 



B. Neutrinos 

Although there are only upper limits determined for 
Era,/ there is clear evidence from neutrino oscillations 
that the neutrinos have mass. Therefore it is more mo- 
tivated to include f v or £m„ in the fits than any of the 
parameters discussed in the previous subsection. There 
is neutrino information in all papers except P5 and P2; 
the latter notes that n s has the same effect as massive 
neutrinos. P3 and P4 obtain relatively small values for 
n s and higher neutrino masses, but this is not a very 
significant correlation. 

The Ly-a forest data present more difficulties than the 
other data sets in the form of nuisance parameters which 
have to be marginalized over and systematic errors which 
have to be estimated. The Ly-a forest data are orthogo- 
nal to the other data in the sense that they are responsi- 
ble for much of the improvement on the primordial spec- 
trum shape and amplitude: as is more accurately deter- 
mined and the neutrino mass limit is much tighter. PI 
admits that more work is needed, reflecting the sceptical 
attitude in the other papers. 

The 95% C.L. results for Sm„ in units of eV are 
PI: < 0.42 (< 1.54 without Ly-a), P3: < 1.07, P4: 
< 1.16. 

Our recommendation is cautiously 



Sm„ < 1.1 eV 



(4) 



In the following we shall use those parameter values 
in PI, P3 and P4 that were obtained with Em^ as the 
seventh free parameter. The P2 data, obtained with 
Em„ = 0, shall be used in a special way to be described 
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shortly. One notes in the PI, P3 and P4 fits that some 
of the parameters are quite insensitive to whether Em^ 
is free or not: this is true for u>b, w m , n s and A s . For 
these parameters we consider the values obtained in P5 
to be sufficiently unbiased to be used here. 



C. Wj,, fib, OJmi fim, h 

Let us compare PI Q and P2 2]. The SDSS "Main 
sample" in P2 is not exactly the same as what PI calls 
" all" . The effects of the spectroscopic LRG data on the 
Main sample are quoted, generally they are small. For 
flat w — —1 cosmology the LRG data change the param- 
eters obtained from the Main sample by the amounts 
Auj m = -0.004, AO m = -0.007, Ah = -0.004, An s = 
-0.017. 

Assuming that the LRG data would cause the same cor- 
rections to the PI results, we concoct a set of LRG- 
corrected results that we call Pl+LRG. 

The LUb values found in the five analyses are very con- 
sistent, uncorrelated with other parameters and robust 
against different choices of parameter sets and priors, as 
pointed out by P4 They are 

PI: 2.36 ± 0.09, P3: 2.24±g'g|, P4: 2.24±^H, P5: 2.31 ± 
0.13. 

We can combine these with the BBN value [2l] 2.2 ± 
0.2, to obtain values of the median and errors 



10 2 w 6 = 2.28±g:i| 



(5) 



The parameter oj m is determined better than f2 m = 
u) m /h 2 , thus the route to fl m is to use u> m and h. The 
values for uj m have often been obtained in this way. Also, 
since WMAP constrains u> m rather than £l m , one then 
avoids a large (Q TO , h 2 ) correlation. We turn first to the 
Hubble parameter h which is measured to be 



Pl+LRG: 0.706 



+0.023 
-0.022' 



P3: 0.648 



+0.039 
-0.038 



, P4: 0.691±0.038. 



P3 and P4 exhibit clearly the effect of making the neu- 
trino mass variable, a change of —3.1% and —6.4% re- 
spectively. This is a measure of neutrino mass bias in 
P2 and P5 that we therefore do not include. The chosen 
input yields the median and errors 



^ = 0.6871°;°$ 



(6) 



in flat space, in remarkably good agreement with our 
constrained fit in 1998 h = 0.68 ± 0.05, based on 
Cepheid distances only. 

From the above values of u>b and h we derive the frac- 
tion of baryonic energy in the Universe 



Q b = 0.048; 



-0.005 
-0.004 



(7) 



Turning now to the m and dm parameters, the avail- 
able input data are the Pl+LRG value f2 m = 0.277 ± 
0.025, the P2 value u m = uj dm + uj b = 0.142 ± 0.005, and 
the Udm values 

P3: 0.126 + 0.007, P4: 0.110 + 0.006, P5: 0.112 + 0.011. 



Combining these values with h and u>b in the most 
efficient way, we find the median and errors 



= 0.139 + 0.011, n m = 0.286 



+0.030 
-0.028- 



(8) 



This is in excellent agreement with what WMAP quoted 
in 2003, Q m = 0.27 ± 0.04 @. It follows from this that 
the dark energy content in flat space is 



= 0.714±°; 028 



030 ' 



and the fraction of dark matter 



fiV/m — 0.238 



+0.030 
-0.028 



(9) 
(10) 



D. 



A s , as 



The parameters r and f2fc are significantly correlated 
so if fife is taken to be zero, also r will be small. PI 
and P2 use the prior r < 0.3 which appears a bit tight; 
even so PI obtains the highest value for r. Since all pa- 
pers use WMAP data, they inherit the problems with r 
appearing different in the Northern and Southern hemi- 
spheres [3 H3] • Thus this parameter is neither reliably 
nor precisely determined. The values that can be quoted 
are 



PI: 0.185 



+0.052 



P3: 0.108 



-0.049 



P4: 



0.1431°;°™ P5: 



-0.046' x "- iuu -0.047' 

0.147+ 0.085. twocolumn From this, the values of the 
median and errors are 



0.147; 



-0.068 
-0.064 



(11) 



The parameter n s is measured by PI, P3, P4 and P5. 
The LRG data in P2 do not yield independent informa- 
tion on n s and u> m , so P2 adds them in the form of a 
constraint uo m = 0.130(n s /0.98)~ 12 ± 0.011 that we do 
not use. The available data are then 
Pl+LRG: 0.972l° ° 2 ^, P3: 0.95 ± 0.02, P4: 0.957^ ^9, 
P5: 0.975 + 0.038.' 

From this, the values of the median and errors are 



n s = 0.962l°;°|? , 



(12) 

thus n s is marginally less than 1.0. 

There are three measurements of the amplitude A s , 
quoted in the form ln(10 10 A s ): 
P3: 3.1 + 0.1, P4: 3.11^;^, P5: 3.17 + 0.16. 



From this, the values of the median and errors are 



H10 W A S ) = 3.12±g;i| 



(13) 



The final parameter in this survey is as- P2 has fixed 
its value at 0.85, the other experiments quote values ex- 
hibiting a rather large variance, 



+0.035 
-0.033' 



P3: 0.74 ± 0.08, P4: 0.6781°,;°^, P5 



PI: 0.890 
0.849 + 0.062. 

From this, the values of the median and errors are 



rr — nS1+ 009 

(Is — U.81_ 14 



(14) 



Since ag is rather strongly correlated with r and fl m , 
a more accurate determination of Q m could have been 
obtained if all experiments had agreed to fix ag and r at 
some common values. 
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TABLE II: Recommended values of parameters 



Parameter 


Definition 


Median value 


83-/17-percentile, 
or 95% CL limit 


10 2 CJ6 


10 2 x baryon density 


2.28 a 


+0.12/ - 0.13 


t~\ / 7 2 

J2(, = LJb/h- 


Normalized baryon density 


0.048 c 


+0.005/ - 0.004 




Total matter density 


0.139 " 


±0.011 


Q>m = OJm/h 2 


Normalized matter density 


0.286 c 


+0.030/ - 0.028 


Qdm — 


Normalized dark matter density 


0.238 c 


+0.030/ - 0.028 


Qa = 1 — — 


Normalized dark energy density 


0.714 c 


+0.028/ - 0.030 


h 


Hubble parameter [100 km/s Mpc] 


0.687 a 


+0.034/ - 0.047 




Neutrino mass sum 




< 1.1 eV 


T 


Thomson scattering optical depth to decoupling 


0.147 a 


+0.068/ - 0.064 


n s 


Scalar spectral index 


0.962 a 


+0.030/ - 0.027 


ln(10 10 A,) 


Scalar fluctuation amplitude 


3.12 a 


+0.14/ - 0.12 


0"8 


RMS linear mass perturbation in 8/t -1 Mpc spheres 


0.81 c 


+0.09/ - 0.14 


r = At/ A s 


Ratio of tensor to scalar amplitude fluctuations 




< 0.4 b 


a s = dn s /dk 


Running scalar index 


-0.011 b 


±0.012 


((! 


Dark energy EOS 


-0.92 b 


+0.17/ - 0.12 




Normalized vacuum density 


-0.023 b 


+0.017/ - 0.050 



a From 7-parameter fits including Em„. 

b From alternative 7-parameter fits which exclude Em, 

c Alternative or derived parameter. 



VI. CONCLUSIONS 

By combining results from five large data analyses 0, 
S S LI we find that no significant values can be obtained 
at present for r, a s , w, Wq, w\, and that the values for 
flk and n s are only marginally significant, consistent with 
flat space and no tilt. This agrees well with a Bayesian 
study [2^ (using much of the PI data) which concludes 
that w and n s do not improve any fits. 



We consider it important to include Hm u among the 
varied parameters since neutrino oscillation results have 
proven that neutrinos have mass. Yet the fits we quote 
yield only an upper limit for Em„. 

We have shown that it is possible to recommend best 
values for u>b, ui m , h, r, n s , A s , erg, and n s , rather 
regardless of whether r and as are taken to be fixed or 
not; fixing would improve the accuracy of several param- 
eters. Our results are collected in Table II. 
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